A first level 
terdisciplinary 


course 


ay 
eal 
taal 
Qa 
Oo 
= 
VW) 
— 
O 
— 
ra 
few 
2 
O 
VU 


ae, Te ac ag ene oe eer iat ee 


oe 


| ae |, Mees, ar : : oa eee 
A : ; . : re . : Cae ' a id > 7 Y 
% is £2) sss 5 en ee ae ek eee 


| BLOCK C 
CONTINUOUS MODELS 


atl 


Mm 


=A 
= 


sing 
c 


© 


Prepared by the course team 


About this course 


This computer book forms part of the course MST121 Using Mathematics. 
This course and the courses MU120 Open Mathematics and MS221 Exploring 
Mathematics provide a flexible means of entry to university-level mathematics. 
Further details may be obtained from the address below. 


MST121 uses the software program Mathcad (MathSoft, Inc.) and other 
software to investigate mathematical and statistical concepts and as a tool in 
problem solving. ‘This software is provided as part of the course. 


This publication forms part of an Open University course. Details of this and 
other Open University courses can be obtained from the Course Information and 
Advice Centre, PO Box 724, The Open University, Milton Keynes, MK7 6ZS, 
United Kingdom: tel. +44 (0)1908 653231, e-mail general-enquiries@open.ac.uk 


Alternatively, you may visit the Open University website at 
http: //www.open.ac.uk where you can learn more about the wide range of 
courses and packs offered at all levels by The Open University. 


To purchase a selection of Open University course materials, visit the webshop 
at www.ouw.co.uk, or contact Open University Worldwide, Michael Young 
Building, Walton Hall, Milton Keynes, MK7 6AA, United Kingdom, for a 
brochure: tel. +44 (0)1908 858785, fax +44 (0)1908 858787, e-mail 
ouwenqg@open.ac.uk 


The Open University, Walton Hall, Milton Keynes, MK7 6AA. 
First published 1997. Second edition 2003. 
Copyright © 2003 The Open University 


All rights reserved; no part of this publication may be reproduced, stored in a retrieval system, 
transmitted or utilised in any form or by any means, electronic, mechanical, photocopying, 
recording or otherwise, without written permission from the publisher or a licence from the 
Copyright Licensing Agency Ltd. Details of such licences (for reprographic reproduction) may be 
obtained from the Copyright Licensing Agency Ltd, 90 Tottenham Court Road, London WIT 4LP. 


Edited, designed and typeset by The Open University, using the Open University TRX System. 
Printed and bound in the United Kingdom by the Burlington Press (Cambridge) Limited. 
ISBN 0 7492 5558 7 

rae 


Guidance notes 


Chapter Cl | 
Section 5 Optimisation with the computer 
5.1 Finding derivatives i 
9.2 Optimisation 
Chapter C2 
Section 5 Integration with the computer 
5.1 Finding indefinite integrals 
5.2 Definite integrals, areas and summations 


Chapter C3 
Section 4 Differential equations with the computer 


4.1 Direction fields and solution curves 
4.2 Euler’s method 


Solutions to Activities 


Index 


_ Guidance notes 


This computer book contains those sections of the chapters in Block C 
which require you to use Mathcad. Each of these chapters contains 
instructions as to when you should first refer to particular material in this 
computer book, so you are advised not to work on the activities here until 
you have reached the appropriate points in the chapters. 


In order to use this computer book, you will need the following Mathcad 
files. 


Chapter Cl 
121C1-01 Differentiation 


Chapter C2 
121C2-01 Integration as a limit of summations (Optional) 


Chapter C3 


121C3-01 Direction fields and solution curves 
121C3-02 Direction fields, solution curves and Euler’s method 


Instructions for installing these files onto your computer’s hard disk, and 
for opening them, are given in Chapter AO. 


The computer activities for Chapters Cl and C2 require you also to work 
with Mathcad documents which you have created yourself. 


Activities based on software vary both in nature and in length. Sometimes 
the instructions for an activity appear only in the computer book; in other 
cases, instructions are given in the computer book and on screen. 
Feedback on an activity is sometimes provided on screen and is sometimes 
given in the computer book. 


For advice on how each computer session fits into the suggested study 
patterns, refer to the Study guides in the chapters. 


Chapter C1, Section 5 
Optimisation with the computer 


The Mathcad activities for Chapter C1 have only one prepared file 
associated with them. This file explains how to differentiate with the 
computer, which is the topic of Subsection 5.1. To solve optimisation 
problems, in Subsection 5.2, you will be given guidance on how to create 
your own Mathcad documents. 


By the end of this section, you should be able to use Mathcad to carry out 
the process of differentiation and to solve optimisation problems, and also 
be a more skilled and independent Mathcad user. 


5.1 Finding derivatives 


In Mathcad, derivatives are found using the < operator. ‘This operator 
can be used in two different ways, depending on whether you simply wish 
to find the numerical value of the derivative of a function f at a particular 
point, or whether you seek a general formula for the derivative of f. 


For example, the derivative of the function f(a) = :x* can be described by 

the formula f’(a) = 2a, which holds for all real values of . Mathcad can 

obtain this expression for the derivative, using the - operator symbolically, 

as indicated in Activity 5.2. The value of the derivative f’(x) = 2x at the 

particular point x = 3 can then be found as f’(3) = = x 3 = 1.2. However, 

it is also possible to obtain the particular value f’(3) = 1.2 without first Essentially, this involves a 
finding the symbolic expression 2x for the derivative. Mathcad can achieve direct application of the 


this by using the operator numerically, as you will see in Activity 5.1. definition of derivative, as 
‘ given by equation (1.3) in 


Chapter C1 (with, in this 


Activity 5.1 Evaluating a derivative numerically a Fad ti = 1 
Open Mathcad file 121C1-01. Page 1 introduces the document. Work Don’t forget to make your 
through page 2, and then carry out Task 1 on page 3. own working copy of the file. 


Solutions are given on page 31. 


Comment 


© It is a good idea to adopt a systematic approach when finding the two 
points at which the derivative of f(x) = x° — x? — 3x + 4 is zero. For 
example, try a value that looks reasonable from the graph, say 
x = 1.5, which gives the derivative value 0.75. This value is positive 
(as is also apparent visually from the slope of the corresponding 
tangent plotted in blue on the graph). Since the value of x sought here 
is a minimum of the graph of f, the positive derivative indicates that 
x = 1.5 is too large an estimate. Hence try further values for z, each 
0.1 less than the previous one, until (for x = 1.3) the value given for 
the derivative is negative. This indicates that a value for x at which 
f'(x) = 0 is somewhere between the last two zx-values tried, 1.3 
and 1.4. If a greater degree of accuracy is required, then this process 
can be repeated by trying further values of x between these last two, 
and continuing until the desired level of accuracy is obtained. 


© 


The main aim of this activity was to illustrate the Mathcad facility for 
evaluating derivatives numerically at given points, though it should 
also have reinforced the link between the value of the derivative at 


Remember that Mathcad 
notes are optional. 


You should still be working 
with Mathcad file 121C1-01. 


These commands are 
introduced in Chapter AO, 
file 121A0-05. 


CHAPTER C1 


various points and features of the graph of the function. In order to 
find where the function f(x) = x° — x* — 3x + 4 has derivative zero, it 
is, in practice, quicker to find the formula for the derivative f’(x), set 
it equal to zero, and solve the resulting equation, 3x7 — 2x — 3 = 0. 


Mathcad notes 


© When evaluating a derivative in this way (numerically), you must 
define earlier in the document the point at which the derivative is to 
be found, e.g. x:= 3. You can then either enter an expression directly 
into the derivative operator, e.g. 27/5, or enter the name of a function, 
e.g. f(x), which has been defined previously. (These two approaches 
are shown on pages 2 and 3 of the document, respectively. ) 


© When you type the equals sign = to evaluate an expression 
numerically, it doesn’t matter where on the expression the blue bar 
cursor is, or if all or part of the expression is enclosed by the blue 
selection box. All that matters is that the expression is complete, with 
every placeholder filled in. 


© When an expression involving the - operator is evaluated numerically, 
Mathcad uses a numerical method to obtain an approximation to the 
exact value of the derivative at a point, which is usually accurate to 7 
or 8 significant figures. Very occasionally the method fails, in which 
case Mathcad displays the error message ‘not converging’. 


Once the symbolic processor is loaded, you can use Mathcad to 
differentiate a function symbolically and hence obtain an algebraic 
expression for the derivative. This replicates what you might do by hand, 
but Mathcad can be used to differentiate functions that would be rather 
complicated to do by hand. 


Activity 5.2 Differentiating symbolically 


Work through page 4 of the document to load the symbolic processor and 
set up how the results are displayed. Then follow the instructions to 
differentiate the function x?/5 symbolically. 


Carry out Task 2 on page 5 of the document. 


Solutions are given on page 31. 


Comment 


© Each of parts (b)-(h) in Task 2 concerns a function that you were 
asked to differentiate in the main text. In the solutions, where the 
expression obtained using Mathcad differs in appearance from that 
obtained by hand, both expressions are given. If you cannot see 
immediately why such a pair of expressions are equivalent, it is a good 
idea to copy down one expression and use algebra to verify that it can 
be rearranged to give the other. 


© Sometimes the form of solution obtained by Mathcad can be 
‘improved’ by simplifying it, as follows. After differentiating 
symbolically, leave the derivative selected (enclosed in the blue 
selection box) and choose Simplify from the Symbolic menu. In 
many cases Factor Expression, also from the Symbolic menu, can 
be used to simplify an expression. However, these options should be 
used with caution, as what Mathcad regards as ‘simpler’ may not 
necessarily seem so to a human observer! 
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: SECTION 5 Optimisation with the computer 


Make sure that the variables entered in the two placeholders of the 
derivative operator match. For example, if you mistakenly try to 
evaluate the expression 


d 
—cos(4- 2), 


dt 


then you will obtain the answer 0. 


Brackets can be crucial when entering an expression into the 
right-hand placeholder of the derivative operator. For example, if the 
expression x° — x” — 32 + 4 is entered without being enclosed in 
brackets, Mathcad will differentiate the x* term but not the others. 


You may find that Mathcad will not allow you to enter into the 
placeholder a quotient whose numerator is a sum of terms (it depends 
on how you enter the quotient). This difficulty can be avoided by 
enclosing the quotient in brackets. 


To avoid such problems, it is good practice always to include outer 
brackets around expressions that you enter into the right-hand 
placeholder of the = operator. This also helps to make the expression 
to be differentiated clear on the screen. 


Mathcad notes 


o, 


When Mathcad carries out symbolic calculations, it does not refer to 
any variable or function definitions that occur in the document. For 
example, symbolic evaluation of the expression 


d 
ay 6 -x°) yields 2-a-z, 
x 


whether or not a has been defined. Similarly, symbolic evaluation of 
the expression 


d d 
whether or not the function f has been defined. In this case, Mathcad 
is unable to evaluate the expression and simply reproduces it. 


Both the above expressions can be evaluated numerically in the usual 
way, provided that the variables x and a, and the function f, have 
been defined previously in the document. 


The expressions e” and exp(zx) are equivalent in Mathcad, but the 
latter form is always used in the output of symbolic calculations, 
irrespective of the form used for input. 


Mathcad carries out a symbolic calculation only when you select an 
expression and choose a command from the Symbolic menu or use 
the corresponding keyboard shortcut. The result of a symbolic 
calculation is not updated if you change the original expression (unlike 
the result of a numerical calculation). 


You may have noticed the Mathcad command Differentiate on 
Variable in the Symbolic menu. This provides an alternative 
approach to finding a formula for the derivative of a function, which is 
described in A Guide to Mathcad. However, it is usually better to use 
the + operator and Evaluate Symbolically, as described in the 
current activity, because this results in clearer Mathcad documents. 


Now close file 121C1-01. 


An example of such a 
quotient is 


sin(t) — t? 
et : 


See Chapter C1, 
Subsection 2.3. 


See Computer Book A, 
Chapter A3, Subsection 5.5. 
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5.2 Optimisation 


Optimisation involves finding the greatest or least value taken by a 
function on an interval. In the main text you saw how to apply the 
Optimisation Procedure by hand, but all the calculations required can also 
be performed using Mathcad. In some cases, Mathcad can be applied to 
check calculations already done by hand, while in other cases, calculations 
can be carried out that are too complicated to do by hand. 


In this subsection you will solve two optimisation problems, the first a 
minimisation problem and the second a maximisation problem. 


The orienteer’s problem 


The orienteer’s problem was described in Computer Book A. The orienteer 
starts at a point in a forest and needs to reach a final point on a path. The 
problem is to find where the orienteer should aim to join the path, where 
the joining point is « km from a fixed point O on the path, in order to 
minimise the time taken overall. For the particular data given, this 
involves finding the value of x that gives the least value of the journey time 
(in hours) 


f(x) =0.125V1 + x? + 0.0625(2— x) (a in [0,2]). 
The graph of y = f(x) is shown in Figure 5.1. 


Figure 5.1 Graph of the function for the orienteer’s problem 


In Computer Book A, you found an approximate solution to this problem 
using Mathcad, by zooming in on the graph and using the crosshair. You 
are now in a position to solve the problem more accurately, using 
differentiation. 


In Activity 5.3 you are asked to use the Optimisation Procedure to solve 
the orienteer’s problem. In place of a prepared Mathcad file for this 
activity, you are guided through creating your own Mathcad document. 
This includes a requirement to enter some explanatory text, to make the 
document more comprehensible. 


The instructions for this activity may look lengthy, but many of them 
describe general Mathcad techniques that can save you time and effort, 
now and in the future. Most of the instructions are given initially in a 
mouse/click way, while details of keyboard alternatives are provided in the 
first Mathcad note for the activity. 


SECTION 5 Optimisation with the computer 


Activity 5.3 Minimising the orienteer’s journey time 


The following instructions guide you through creating a Mathcad 
document to carry out the Optimisation Procedure, in order to solve the 
orienteer’s problem. Part (b) describes how to enter a title in the 
document. You should similarly enter any other text that you think is 
appropriate. For example, you could include a line of text before each step 
of the Optimisation Procedure, to explain what you are doing. 


(a) Begin by creating a new blank document, as follows. With Mathcad 
running on your computer, select the File menu and choose New. (If 
you have just started Mathcad running, then there is no need to do 
this, as it automatically starts with a new document.) 


Enter a title at the top of your document. To do this, click to position 
the red cross cursor in an appropriate place, and choose Create Text 
Region from the Text menu. Then type a suitable title, for example, 
Optimisation - the orienteer’s problem, in the text box. To 
finish, click anywhere outside the text box. If you need to edit the text 
later, simply click on it. 


Load the symbolic processor by choosing Load symbolic processor 
from the Symbolic menu. (You will be using both symbolic 
differentiation and symbolic solution of an equation in the 
optimisation process.) | 


It is also helpful to use the Mathcad facility for providing comments, 
by selecting Derivation Format... from the Symbolic menu, 
clicking in the check box to ‘Show derivation comments’ and selecting 
the option button to show the derivation steps ‘vertically, inserting 
lines’, then clicking ‘OK’. 

Now you are ready to solve the orienteer’s problem, by carrying out 
the three steps of the Optimisation Procedure, as follows. 


Step 1: Find the stationary points of f 


Use the £ icon on palette 1 and Evaluate Symbolically from the 
Symbolic menu to differentiate the expression to be optimised, 
which is 


0.125-V1+ 2? + 0.0625 - (2 — 2). 
You can use the key sequence 
(0.125*\1+xA2([t] [1] [1] [7] +0.0625* (2-x) ) 


to enter this. The \ (backslash) keystroke creates the square root sign. 
Remember to enclose the expression to be differentiated in outer 
brackets. 


To find the stationary points, select an occurrence of x in the 
expression you have just found for the derivative (put the blue bar 
cursor next to one of the zs) and choose Solve for Variable from the 
Symbolic menu. You should see two solutions, —0.577 350... and 
0.577 350..., within a column vector. These are the values of x where 
the derivative is zero; that is, the stationary points. 


A Guide to Mathcad contains 
detailed information on 
creating and editing your own 
documents. 


Any text that you enter can 
later be moved, or deleted, as 
you wish. 


If you have already loaded the 
symbolic processor during this 
Mathcad session, then you 
will not need to do so now. 


The derivation steps were 
shown horizontally in 
Activity 5.2, but here they 
should be shown vertically. 


If you have difficulty in 
following these instructions, 
then you may like to look 
ahead to the first note in the 
Comment below, which shows 
what should eventually 
appear on the screen. 


The derivative should appear 
underneath the original 
expression. If not, then 
review the instructions above 
about ‘Derivation Format...’. 


To obtain the blue bar cursor, 
click on the expression 
(anywhere will do), then press 
the ‘Down Arrow’ key [J] 
repeatedly until the cursor 
appears. You can use the left 
[+] and right [—] arrow 
keys to move the cursor into 
position next to an z. 


ee 


Step 2: Evaluate f at each of the relevant points 


Use Mathcad to evaluate the original expression at the two endpoints 
of the interval, zc = 0 and x = 2, and at the stationary point between 
them. 


There are several ways of doing this. Probably the most efficient is to 
define the function 


f(x) = 0.125 - V1 4+ x? + 0.0625 - (2 — 2). 
You can then evaluate the function for any given value of z. 


Rather than create the whole function definition from scratch, you can 


Copy and paste is a handy take advantage of Mathcad’s copy and paste facilities. You need to 
way of avoiding the need to create the left-hand side of the definition, f() := (that is, type 

retype awkward or lengthy f(x):), but then, for the right-hand side, you can copy the expression 
expressions, so it is worth that you differentiated in Step 1. To do this, select it (click anywhere 
remembering. in the expression and press the ‘Up Arrow’ key [*] until the whole 


expression is enclosed within the blue selection box) and choose Copy 
from the Edit menu. Then paste this expression into the function 
definition, by first clicking on the placeholder following f(z) := to 
select it, and then choosing Paste from the Edit menu. 


Once you have defined the function f(x), evaluate f at the interval 
endpoints, z = 0 and x = 2, and at the stationary point; for example, 
type f(0)= to evaluate f(0). It is sufficient to input the value for the 
stationary point to three decimal places, x = 0.577, or you can use 
copy and paste to input all twenty digits given in the solution if you 
wish! 


Step 3: Choose the optimum value 


Remember that the Identify the least of the three function values that you found. This is 
orienteer’s problem is a the minimum value of f(x) within the given interval. 


minimisation problem. Bie 
You could finish at this point, but it is a good idea to record your 


conclusions in the document. For example, you could enter the text 
The least value of f(x) on [0,2] is ? (at x=7). 


replacing the question marks with the values that you found. You can 
add more text if you wish, but the Optimisation Procedure for the 
orienteer’s problem is now completed. 


(e) Finally, save your file, by choosing Save As... from the File menu. 
You will need to give it a suitable name that indicates the contents, for 
example, my121C1-orienteer. This will help you to organise and 
identify your files. 


Solutions are given on page 31. 


+ SECTION § Optimisation with the computer 


Comment 


© The Mathcad document should now look something like this: Your document should 
contain the calculations 


shown but may have different 
Find the distance x that minimises the time taken to run f(x). text. 


Optimisation - the orienteer's problem 
Step 1 Find the stationary points of f. 

d 2 

ee 14x" +0.0625-(2 - x) 
dx 

yields 


eg ee 


2 
+z 


has solution(s) 


- .37735026918962576451 
37735026918962576451 


Step 2 Evaluate f at the endpoints x = 0 and x = 2, and at the stationary point x =0.57? 


f(x) :=0.125-9[1 + x7 + 0.0625-(2- x) 


f(0) = 0.25 f(2) = 0.28 (0.577) = 0.233 


Step 3 Choose the optimum value 
The least value of f(x) on [0,2] is 0.233 (at x = 0.577). 


© When an expression does not contain an equals sign, Solve for 
Variable finds the values of the selected variable that make the 
expression equal to zero. Note that Mathcad omits the zero before the 
decimal point in the solutions. 


Mathcad notes 


© Mathcad provides keyboard alternatives for most of the menu 
commands and features used above, as follows: 


new blank document type [F7] (function key ‘F7’, not the 
separate keys ‘F’ and ‘7’); 
“ icon (on palette 1) type ? (a question mark, given by [Shift] /); 


create text region type " (a double-quote, given by [Shift] 2); 
evaluate symbolically type [Shift] [F9]; 

copy type [Ctrl]c; 

paste type [Ctrl] v. 


© When working symbolically, a decimal point in the input expression 
triggers a decimal result, with up to twenty decimal places! However, 
when working numerically (evaluating an expression using =), the 
number of decimal places displayed depends on the ‘Displayed 
Precision’ (Math menu, Numerical Format...) specified for the 
result. (By default, the results of numerical calculations are displayed 
to three decimal places.) 


Now close the Mathcad file that you have created and saved. 


J GHAPTERTI 2 


A traffic planning model 


Traffic planners wish to set up a mathematical model to describe how the 
volume flow rate of traffic (that is, the number of vehicles which pass a 
fixed point in a given time) varies with the average velocity of vehicles 
along a single lane of a road. The eventual purpose of this model is to 
advise on how traffic flow can be maximised. 


The planners assume that each vehicle moves at a constant velocity 

v ms‘. On the basis of a subsidiary model and many observations, they 
estimate that, on average, each driver maintains a distance v + 0.02v’ 
metres between the front of their vehicle and the back of the vehicle 
immediately ahead. The average length of a vehicle is estimated to be 5 
metres. ‘The model therefore represents the situation as shown in 


Figure 5.2. 
velocity v 
ie J 
v + 0.02v? v + 0.02v? 


5 is) 5 
Figure 5.2 Representation of traffic flow 
The distance between the fronts of successive vehicles is 5 + v + 0.02v” 
metres. Each vehicle, travelling at velocity v ms~', covers this distance in 


(5+ v+0.02v?)/v seconds. It follows that the volume flow rate f(v) of 
traffic (in vehicles per second past a fixed point) is given by 


Vv 
ee eee 
A velocity of 35 ms~? is This formula is to apply for 0 < v < 35. 


about 126 kilometres per hour 


The traffic planners seek the maximum value of f(v) for v in the interval 
or 78 miles per hour. 


(0, 35]. If such a maximum can be found, it can be used to provide an 
advised speed of travel on the road. 


In the next activity you are asked to create a Mathcad document once 
more, to find both the greatest value of the traffic flow rate function f(v) 
and the value of v for which this occurs. 


Activity 5.4 Maximising the traffic flow rate 


You do not need to load the (a) Create a new blank document, load the symbolic processor, and set 
symbolic processor or set the the derivation format as described in Activity 5.3(c). 
oe ee ee (b) Enter a title at the top of your document, for example, Maximising 


continuing directly from 


Activity 5.3. the traffic flow rate. 


(c) Solve the traffic flow rate problem by finding the maximum value of 
the function 


v 
v) = ——___~ 
P(e) 5+v+0.02v? 
Remember that this is a on the interval [0,35] and the corresponding value of v. Do this by 
maximisation problem. carrying out the three steps of the Optimisation Procedure, just as you 


did in Activity 5.3(d). | 
(d) State your answer in a form which is appropriate in the context of the 
traffic planning model. 


Solutions are given on page 31. 


aS 


Comment 


If you have defined a function f(v) for the volume flow rate of traffic in 
your document then you may wish to display a graph of it. You can do so 
by defining a suitable graph range, for example, v:= 0,0.1..35, and then 
plotting f(v) against v. While the graph confirms the maximum value at 
v ~ 16, it also shows that there is little difference in the value of f(v) 
between v = 10 and v = 20. 


Save your file, for example as my121C1-traffic. Then close this file. 
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Chapter C2, Section 5 
Integration with the computer 


If necessary, see Chapter C1, 
Activity 5.3(a) and (c) on 
page 9 of this computer book. 
However, note that the 
derivation steps here are to 
be shown horizontally. 


Be careful not to confuse the 
palette 3 icon { with the 


palette 1 icon . which is 
used for finding definite 
integrals (as you will see 
later). 


It is a good idea to insert a 
title at the top of any 
document that you save. If 
you need to create space for 
this, then do so by putting 
the red cross cursor at the top 
of the document and using 
Ins/Del Blank Lines... 
from the Edit menu. 
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There is only one prepared Mathcad file for this section, and that comes 
towards the end and is optional. As in the later computer activities for 
Chapter C1, you will, for the most part, be creating your own documents 
here and using Mathcad directly. Subsection 5.1 shows how to find 
indefinite integrals in Mathcad, while Subsection 5.2 covers definite 
integrals. 


5.1 Finding indefinite integrals 


In Mathcad, indefinite integrals are found using the { operator. Unlike the 
= operator, the { operator can be used only symbolically, to find an 
algebraic expression for an integral of a given function. It cannot be used 
numerically. 


In this subsection you are invited to find indefinite integrals for a variety of 
functions. The first activity provides an introduction to finding indefinite 
integrals using Mathcad. 


Activity 5.1 Indefinite integrals 


In this activity you will use Mathcad to find the indefinite integral of x’. 


(a) Create a new blank Mathcad document and load the symbolic 
processor. Then set the derivation format to show derivation 
comments and to display the derivation steps horizontally. 


(b) Enter the { operator in your document, either by clicking on the 
palette 3 icon {, or by using the keyboard shortcut [Ctr1l]i. (To 
reveal palette 3, click on the number button ‘1’ at the top of palette 1 
on the left-hand side of the Mathcad window, then on the number 
button ‘2’.) 

(c) Enter the expression to be integrated, x”, in the first placeholder after 


the integral sign. (This expression is called the integrand.) Then enter 
the variable of integration, x, in the placeholder after the ‘d’. 


(d) Use the ‘Up Arrow’ key [1] to select the entire expression (including 
the integral sign). Then choose Evaluate Symbolically from the 
Symbolic menu, or use the keyboard shortcut [Shift] [F9]. Check 
that the answer provided by Mathcad is what you expect. 


(e) Now go through the same procedure to evaluate the integral / u? du. 


(f) If you wish to save your work, then select the File menu and use Save 
As... to name and save your document. 


SECTION 5 Integration with the computer 


Comment 


© Notice that the f operator in Mathcad gives only an integral of the 
integrand. It does not give the indefinite integral because it does not 
add an arbitrary constant. 


© The outcomes from integrating fx? dz and f{ u? du demonstrate that 
the form of the indefinite integral depends on the nature of the 
function being integrated but not on the choice of symbol for the 
variable of integration. 


Mathcad notes 


© Remember that symbolic evaluation occurs only when you select 
something and choose the command from the menu, or use the 
keystroke [Shift] [F9]. The result of a symbolic evaluation is not 
updated automatically if you just change the original expression. In 
part (e), for example, it is not sufficient simply to edit the integral 
expression { x? dx, replacing each ‘2’ by a ‘u’. 


© If you enter the expression { x? dx and then press the = key, the error 
message ‘Symbolic operator only’ appears, whether or not you have 
previously specified a value for x. This is as it should be, since the = 
key is used for numerical evaluation, which does not make sense in the 
context of finding an indefinite integral. 


© You cannot find an integral of a function f by first defining f and then 
evaluating the expression [ f(x) dz. If you try it, then Mathcad will 
simply repeat the expression. There has to be an actual algebraic 
expression in the integrand placeholder for symbolic integration to be 
carried out. 


In each of the remaining activities in this subsection, you can either 
continue using the document from the previous activity, or close that file 
and create a new blank document. You will need to decide whether you 
wish to save your work. 


For each activity, you will also need to ensure that Mathcad is prepared to 
carry out symbolic calculations. The symbolic processor should be loaded, 
and the derivation format set to show derivation comments and to display 
the derivation steps horizontally. 


In the next activity, you are asked to use Mathcad to find the indefinite 
integrals of two functions that you integrated by hand in the main text, 
and of a third function which you could also integrate by hand. 


Activity 5.2 Finding indefinite integrals 


Use Mathcad to evaluate each of the following indefinite integrals. (If 
necessary, refer to the instructions for evaluating an integral in 
Activity 5.1(b)—(d).) Enter each integral on’a new line. 


ba IG +e) fe te ot 
(b) IG + 5sin(5y) dy (y>0) 


(c) / (a + cos(az))dx (where a is a non-zero constant) 


Solutions are given on page 32. 


If you arrive at incorrect 
answers, then see the 
Mathcad note below. 


15 


Note that sin? x should be 
input as sin(x)A2. 


16 


OMPTER 


Comment 


© After the ‘+c’ has been added, the Mathcad answers to parts (a) 
and (b) agree with those obtained in the main text, and that for 
part (c) is identical to the answer found by applying Table 1.1. In 
part (a), recall that exp(3z) is an alternative way of writing e°”. 


© Note that, while Mathcad gives answers here that agree with those 
obtained by hand, there is no way of entering the additional 
constraints on the variables that accompany the integrals, for example, 
the condition xz > 0 in part (a). For this reason, it is necessary to bear 
in mind that the symbolic manipulations performed by Mathcad might 
not be valid in all circumstances, and that the output should be 
interpreted with care. 


Mathcad notes 


You need to type 3*x for 3x in part (a), 5*sin(5*y) for 5sin(5y) in 
part (b) and a*x for az in part (c), because all products must be made 
explicit in Mathcad. For example, if you type 3x rather than 3*x then 
Mathcad will assume that you have entered a single variable name, ‘32’, 
rather than the product of the number 3 and the variable zx. 


The next activity contains four integrals which you evaluated by hand in 
Section 2, followed by three further integrals (in parts (e)-(g)) that cannot 
be evaluated by hand simply on the basis of what is given in the chapter. 


_ Activity 5. 3 | Further integr als 


Use Mathcad to evaluate each of the following indefinite integrals. 


(a) [(e-3)@-1)de (b) [=e « 
(c) [ sin? xdz (d) / Sood 


z?2+1 


(e) [we du (f) [ ?n(5z) dx (g) [yee 


Solutions are given on page 32. 


Comment 


© The answers obtained from Mathcad in parts (a)—(d) are equivalent to 
the answers found by hand. However, they are not always given in an 
identical form, and in part (c) some algebra is required to show the 
equivalence of the two forms. 


o 


It is possible to use the eo operator to differentiate symbolically each 
of the integrals obtained in this activity. This leads to an expression 
equivalent to the original integrand in each case, as expected. 
However, the output obtained is not always identical in form to that of 
the original integrand, and again some algebra is sometimes required 
to show the equivalence of the two forms. 


© 


Recall that Simplify from the Symbolic menu can sometimes be 
used to simplify to some extent an expression obtained using 
Evaluate Symbolically. In many cases Factor Expression, also 
from the Symbolic menu, can be applied in a similar way. 


SECTION 5 “Integration with the computer S 


Mathcad notes 


The Mathcad command Integrate on Variable, available from the 
Symbolic menu, provides an alternative approach to finding a formula for 
an integral of a function. This is described in A Guide to Mathcad. 
However, it is usually better to use the { operator and Evaluate 
Symbolically, because this results in clearer Mathcad documents. 


Parts (e)-(g) of Activity 5.3 show that Mathcad can extend your 
‘integration reach’ beyond the types of functions that you have so far 
learnt how to integrate by hand. However, Mathcad cannot integrate every 
function. For example, if you ask Mathcad to evaluate symbolically either 
of the indefinite integrals 


/ edt and / Vsin x dx, 


the response is for Mathcad to repeat the given integral without alteration 
(except for the replacement of e~* by exp(—t*) in the first case). This is 
how Mathcad responds when it cannot find an integral. This need not be 
due to a shortcoming of Mathcad. In the two cases here, no expression 1s 
known for either integral in terms of standard functions. 


5.2 Definite integrals, areas and summations 


You will need Mathcad file 121C2-01 for (optional) Activity 5.8, later in 
this subsection. First, however, you are invited to create your own 
Mathcad documents as in the previous subsection, but now to find definite 
rather than indefinite integrals. 


In Mathcad definite integrals are found using the 5 operator. This 
operator can be used either symbolically or numerically. If the operator is 
used symbolically, then Mathcad finds an algebraic expression for an 
integral, evaluates this expression at the upper and lower limits of 
integration, and subtracts the second value from the first to find the 
answer. This is the same as the usual approach to finding definite integrals 
by hand. If the operator is used numerically, then Mathcad does not find 
an algebraic expression, but instead uses a numerical method to find an 
approximate value for the definite integral. 


In each of Activities 5.4-5.7 below, as in Subsection 5.1, you can either 
continue using the document from the activity before, or close that file and 
create a new blank document. You will need to decide whether you wish to 
save your work. 


For each activity, you will also need to ensure that Mathcad is prepared to 
carry out symbolic calculations. The symbolic processor should be loaded, 
and the derivation format set to show derivation comments and to display 
the derivation steps horizontally. 


Activity 5.4 introduces you to the symbolic use of the i operator. 


ee 


The [Tab] key provides a 
good way of moving around 
the placeholders. 


In this activity you will use Mathcad to evaluate the definite integral 


3 
1 
[ ia. 
> 


(a) Enter the J operator in your document, either by clicking on the icon 
in palette 1, or by using the keyboard shortcut & (the ampersand sign, 


given on the keyboard by [Shift]7). 


(b) Enter the integrand, the variable of integration, and the upper and 
lower limits of integration in the appropriate placeholders. 


(c) Select the entire expression, including the integral sign, and choose 
Evaluate Symbolically from the Symbolic menu, or use the 
keyboard shortcut [Shift] [F9]. You should obtain the answer 


In(3) — In(2). 


(d) Leave the answer selected, and press the = key. You should obtain the 
answer 0.405, which is the value of In3 — In2 to three decimal places. 


Comment 


Evaluating symbolically an expression involving the i operator gives an 
expression which is an ezact answer (unless the original expression 
contains a decimal point; see the second Mathcad note below). In the 
example in this activity the expression is In(3) — In(2). You can display the 
decimal value of such an expression by evaluating it numerically. This is 
done by selecting the expression, and then pressing the = key. (If the 
expression has just been obtained as a result of a symbolic calculation, 
then it will already be selected.) 


Mathcad notes 


© When you evaluate numerically an expression in Mathcad, the number 
of decimal places displayed is determined by the value of ‘Displayed 
Precision’. The default value of this is 3, but you can change it by 
choosing Numerical Format... from the Math menu. 


© Ifa Mathcad expression involving the f operator has a decimal point 
in any constant in the integrand, or in both limits of integration, then 
evaluating the expression symbolically gives a decimal answer with up 
to twenty decimal places. (Such an answer is unaffected by the value 
of ‘Displayed Precision’.) For example, evaluating symbolically the 


Mathcad expression 


*10 
[ =e yields .40546510810816438198. 
2 4 63 


The next activity provides further practice in using Mathcad to evaluate 
definite integrals symbolically. 


ee — - SECTION 5 Integration with the computer | oF — 


Each of parts (a)—(d) below gives a definite integral that you were asked to 
evaluate by hand in the main text. In each case, use Mathcad to evaluate 
the definite integral symbolically to obtain an exact answer, and then 
evaluate this answer numerically, to display it as a decimal value. 


2 m/4 
(a) / e* dt (b) / (cos(5az) + 2sin(52)) dx Remember that 7 can be 
o obtained from palette 1, or by 
typing [Ctr1]p. 


2 6 7 : 
(c) / va du (d) / e’ sint dt 


Solutions are given on page 32. 


In the next activity you will see an example of a definite integral that 
Mathcad is unable to evaluate symbolically. When this happens it is often 
worth attempting to evaluate the integral numerically, and the activity 
shows you how to do this. 


Activity 5.6 Evaluating definite integrals numerically _ 


(a) Use Mathcad to try to evaluate symbolically the definite integral 


1 
/ e* di. At the end of Subsection 5.1 
0 it was pointed out that 
You should obtain the warning message ‘No closed form found for Mathcad cannot evaluate the 
integral’. This means that Mathcad has been unable to calculate an indefinite integral 
algebraic expression for an integral of the integrand, and so it cannot = 
evaluate symbolically the given definite integral. / e a. 


Click on ‘OK’ to remove the message. 


(b) Evaluate the definite integral numerically, as follows. It should still be 
selected; if it is not, then select it by putting the blue bar cursor 
anywhere within it, or by surrounding all or part of the expression 
with the blue selection box. Then type =. 


You should find that the answer 0.808 is displayed. 


Comment 


© In general, to evaluate a definite integral numerically, you should enter 
it in the same way as for symbolic evaluation, then select it and type =. 


© The reason why some definite integrals can be evaluated numerically 
but not symbolically in Mathcad is that symbolic evaluation requires 
Mathcad to find an algebraic expression for an integral, whereas 


numerical evaluation involves the use of a numerical algorithm. The Activity 5.8 indicates a 
answer obtained from this algorithm is an approximation, though possible basis for such an 
usually an accurate one. . algorithm. 


Mathcad notes 


On rare occasions, the numerical method used by Mathcad for evaluating 
definite integrals fails to produce a value. In such a case, the integral is 
marked with the error message ‘not converging’. 


The next activity gives you further practice in using the Ss operator 
numerically. 


ce 


The first two steps here are 
identical to those for symbolic 
integration. These steps 
would also have been required 
for the definite integral in 
Activity 5.6(b), had the 
expression not previously 
been entered. 


For example, the values of all 
the definite integrals in 
Activity 5.5 can be found 
satisfactorily using numerical 
integration. 


The remainder of this 
subsection will not be 
assessed. 


We assume (as in the main 
text) that f(x) takes only 
non-negative values in the 
interval |a, b]. 


See equation (4.1) in 
Chapter C2, Subsection 4.2. 


“CHAPTER C2 


Activity 5.7 Further definite integrals (numerically) = 


1 
(a) Evaluate numerically the definite integral , t* dt, as follows. 
0 


(i) Either click on the palette 1 icon f or type & to insert the definite 
integral operator and its placeholders. 


(ii) Fill in the four placeholders appropriately. 


(iii) The cursor should be on the expression; if not, select some part of 
the integral. Then evaluate the expression numerically, by typing =. 


1 
(b) Evaluate numerically the definite integral / V1— 2? dz. 
at 


Solutions are given on page 32. 


Comment 


These two definite integrals can be evaluated either numerically or 
symbolically. As you can check, symbolic evaluation gives + for (a) and $7 


for (b). 


You have seen that in Mathcad many definite integrals can be evaluated 
either symbolically or numerically. This raises the question of which it is 
more appropriate to invoke in any given situation. If you want an exact 
answer (so you can see where constants such as 7 feature in it, for 
example), or if a general result is required, such as a formula for 


‘ 
/ cos(ax)dx (where a is a non-zero constant), 
0 


then use the symbolic approach. If you simply want a number that is an 
accurate value for the definite integral, then numerical evaluation should 
suffice. If you want both an exact value and a decimal value for the 
answer, then you can first evaluate the definite integral symbolically and 
afterwards press =. 


Integration as a limit of summations 


In Subsection 4.2, you saw that a definite integral could be approximated 
as closely as required by a finite sum. This was demonstrated in particular 
for the case in which the value of the definite integral gives the area 
beneath the graph of a function, and each finite sum represents the overall 
area of a set of rectangles. Each rectangle is based upon a subinterval, and. 
the approximation to the definite integral improves as the number of 
subintervals is increased. 


For the area beneath the graph of the function f(x) between x = a and 
x = b, which is given exactly by the definite integral 


[ fae, 


the approximation based on N subintervals is 


N-1 


N-1 
—_ 
Shf(atih) =h S~ f(a+ih), where h = ——. 


“SECTION 5 Integration with the computer 


The remaining (optional) activity in this section asks you to use the 
prepared Mathcad file 121C2-01 to explore the relationship between these 
approximations to the area under the curve and the definite integral itself. 


- Activity 5.8 Integration as a limit of summations (Optional) — 


(a) Open Mathcad file 121C2-01 and read through the document, which 
consists of a single page. The definite integral being approximated 
here is 


40 te 
¥ 15,/sin (=) dx, 
0 AO 


whose value was sought by a process of successive approximation in 
Subsection 4.2. 


The value of this integral is the area on the graph beneath the red 
curve. The value calculated for A is an estimate for this area, using 
approximation by rectangles based on N subintervals. The value of NV 
is initially set to 4. 


Investigate the effect of increasing the number of subintervals, N. Use 
in turn the following values for NV: 


20, 50, 100, 500, 1000, 5000, 10000. 


Compare the corresponding values obtained for A with those in the 
right-hand column of the table below. 


Number of Sum of areas 
subintervals of rectangles 

4 402.27 

20 452.71 

50 456.41 

100 457.21 

500 457.62 

1 000 457.64 

5 000 457.66 

10 000 457.66 


(b) You may like to use the file for other definite integrals, to investigate 
the behaviour of area estimates A as N is increased. For example, you 
could investigate these estimates for the definite integral 


1 
Sa 
[ eae, 

0 


by first editing the file so that f(x) =e7* and b=1. 


If you have time, look also at the behaviour of the corresponding area 
estimates for the definite integral 


m/4 
/ tan x dz. 
0 


This is Table 4.1 in 
Chapter C2, Subsection 4.2. 


See the third Mathcad note 
for Activity 5.1 and the first 
Mathcad note for 

Chapter C1, Activity 5.2. 


Comment 


7 


You should find that the numerical values obtained for A in part (a) 
match to two decimal places the values given in the right-hand column 
of the table. These area estimates appear to tend to the limit 457.66 
(to two decimal places), which is also the numerical value provided by 
Mathcad for the definite integral. 


The values of the definite integrals suggested in part (b) are 


1 : m/4 
/ ee” dz = 0.806 suc / tana de = 0.347, 
0 0 


each to three decimal places. In each case, the estimates A converge to 
this value as N increases. 


In the first example, the area estimates are greater than the value 
given by the definite integral for the area under the curve, while in the 
second example they are smaller. You can see this illustrated on the 
graph by setting small values of N (N < 20, say). In the first example, 
the tops of the rectangles all lie above the curve, while in the second 
example they all lie below. 


Mathcad notes 


© 


The summation sign is obtained from the palette 1 icon }~""_, or by 
typing [Ctr1]$ (you have to press the three keys[Ctr1l], [Shift] and 
4 together). 


Note that the definite integral is set up in this file with integrand f(z), 
where the function f(x) is defined earlier in the document. This 
presents no problem for numerical evaluation of the definite integral. 
However, it is not possible to perform symbolic evaluation of a definite 
integral in this form. This is in line with previous Mathcad notes 
concerning symbolic evaluation and functions defined earlier in a 
document. 

The rectangles are filled in by drawing zig-zag lines very close together. 


On the screen these lines give the appearance of a solid block of colour, 
but the tiny gaps between the lines may become apparent if printed. 


Now close Mathcad file 121C 2-01. 


Chapter C3, Section 4 


Differential equations with the computer 


In this section you will use the computer to see how the information 
contained in a differential equation can be displayed graphically, and how 
differential equations can be ‘solved’ numerically and graphically, even 
where no formula for the solution can be found. 


There are two prepared Mathcad files associated with this section. The 
first draws direction fields and plots solution curves (graphs of solutions) 
for first-order differential equations, and the second illustrates how a 
numerical method for obtaining the solution to an initial-value problem 
works in practice. 


This numerical method can be applied to any differential equation of the 
form 


dy 
= = f(z, y), 


where f is a known function of two variables. This form includes the types 
of differential equation considered in the main text, but also others. 


Mathcad does not contain any symbolic facilities for solving differential 
equations directly, to find a formula for the solution. However, as you saw 
in Chapter C2, Mathcad can be used to find integrals, and finding integrals 
is the main constituent of the two methods for solving first-order 
differential equations that are introduced in the main text (direct 
integration and separation of variables). Hence Mathcad can be used 
indirectly to help solve some first-order differential equations. Mathcad 
can also be applied to differentiate a solution that you have already found, 
to check that it is indeed a solution to the given differential equation. 


In this section, however, Mathcad will be used numerically and graphically 
rather than symbolically. 


4.1 Direction fields and solution curves 


The differential equation 


dy 

dx 
cannot be solved, using the methods of this chapter, to give an equation 
relating x and y. However, the direction field of this differential equation 
provides enough information to indicate the different types of solution 
curve that occur. In Activity 4.1 you will see this direction field drawn by 
Mathcad, and will be able to ask Mathcad to plot the solution curve 
through any point of your choice. 


As special cases, f may be a 
function of x alone or of y 


alone. 


Differentiation was the 
subject of Chapter Cl. 


os 


The initial condition is 
y(Zo) = Yo- 
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CHAPTER C3 


Open Mathcad file 121C3-01. The document opens with the direction field 
of the differential equation 
dy 


a 


drawn for a grid of points with integer coordinates in the (2, y)-plane, for 
values of x between —5 and 5, and values of y between —2 and 4. For each 
such point (x,y) a line segment is plotted through the point, and the slope 
of this line segment is f(x,y) =x+ y. 


(a) Set S to 1, so that a solution curve is plotted through the point (0, 0). 
Briefly describe the curve obtained, or make a small sketch of it. Now 
change the value of yo, to obtain a solution curve through each of the 
following points in turn: 


(6.1)... (02) 1), (6 5) 


In each case, note down a brief description of the curve or make a 
small sketch of it. 


(b) Can you group the solution curves seen in part (a) into distinct types 
of behaviour? How many different types of behaviour are there? 


For each of the following points, try to predict from looking at the 
direction field which type of behaviour the solution curve through the 
point will exhibit: 


(—3,-1), (1,0), (4,2). 


Then, by making a suitable choice of values for 7) and yo, use 
Mathcad to plot the corresponding solution curve and to confirm your 
prediction. 


Solutions are given on page 33. 


Comment 


The line segments of the direction field give a good indication of where the 
solution curves lie. 


Mathcad notes 


© A definition for a function of two variables is created in the same way 
as that for one variable. For example, to define the function 
f(x,y) =«x+y in Mathcad, you could type f(x,y) : x+y. 


© The calculations used to draw the direction field and solution curve 
are ‘hidden’ off the page of the Mathcad document to the right. The 
area beyond the right-hand margin and edge of a Mathcad page 
(which are marked by solid and dotted vertical lines, respectively) can 
be used in the same way as the rest of the document. You can place 
mathematical expressions, text or graphs there. You do not need to 
look at the calculations in this document, but in general, you can view 
what is in the ‘hidden’ area by using the horizontal scroll bar to move 
to the right. 


In the next activity you are asked to use Mathcad to plot the direction 
fields for two other first-order differential equations. In each case, you can 
try to visualise from the direction field how the solution curves behave, 
before plotting some of these curves. 


ee ‘SECT, TON 4 Differential equations with the computer — 


Activity 4.2 Direction fields and solution curves 


(a) Investigate the direction field and solution curves for the first-order You should still be working 
differential equation with Mathcad file 121C3-01. 
dy 
a eee 8 
dx 
as follows. 


(i) First set the variable S to zero, so that no solution curve is 
plotted. Then enter the right-hand side for this differential equation 
into the definition for f(z, y). 


(ii) Try to predict from the direction field where the solution curves 
lie. What types of behaviour will the solution curves exhibit? 


(iii) Set the variable S to 1 to display a solution curve, and then try 
different values for x) and yo to confirm the predicted behaviour of 
solution curves. Use your observations to describe the behaviour of the 
solution curves for the differential equation. 


(b) Now follow the same procedure as in part (a) to investigate the 
direction field and solution curves for the first-order differential 
equation 

2 a ae 
However, in this case, start by altering the scope of the grid in the 
y-direction, by setting Y1 := —5, Y2:=5 and q := 10. 


Solutions are given on page 33. 


Comment 


© The differential equation in part (a) is of the form dy/dx = f(x), and 
so can be solved in principle by direct integration. However, it turns 
out not to be possible to find an algebraic expression for 


[ceo pie 


The presence of the cosz means that the slopes of the direction field 
repeat at horizontal intervals of 27. (The slopes are also invariant in 
the vertical direction, for any given choice of x, as noted in the 
solution.) The alternate positive and negative slopes indicate that 
solution curves will undulate. The overall increasing trend is not so 
obvious from the direction field, but might be expected because the 
magnitude (steepness) of the slopes appears to be greater where the 
slopes are positive than where they are negative. 


© The differential equation in part (b) cannot be solved by the methods 
of this course. Using Mathcad to plot the direction field provides a 
good approach to visualising the types of solution curve for such a 
differential equation. 


Now close Mathcad file 121C3-01. 


4.2 Euler’s method 


In Subsection 4.1 you saw how a direction field can be used to visualise the 
information provided by a first-order differential equation. You also saw 
solution curves plotted on top of the direction field. It is straightforward to 
plot such curves for a first-order differential equation whose solutions can 
be expressed in terms of a simple algebraic formula. However, even where 

This was the case for each of | no such formula can be found, approximate solution curves for a 

the differential equations in differential equation can still be plotted. These graphs are based on a 

Activity 4.2. sequence of numerical estimates for solution values, which constitute a 
numerical solution to the differential equation. In this subsection you will 
see how such a numerical solution can be obtained. 


The procedure to be described below, for obtaining an approximate 
numerical solution to a first-order differential equation, is known as Euler’s 
method. ‘To see how the solution is built up step by step, it is illuminating 
to consider the corresponding graphical construction. This involves 
forming a connected chain of line segments, each of which has a gradient 
Recall that a direction field given by the slope of the direction field at the left-hand end of the line 
provides a slope value at segment, as shown in Figure 4.1. 
every point within a given 
region, and not just at the 
particular grid points y 
displayed on a graph. 


slope = f(x,y) 


Figure 4.1 Graphical construction for Euler’s method, where dy/dzx = f(z, y) 


You will see this idea explained in greater detail in the next activity. 


Activity 4.3 Introducing Euler’s method 
Open Mathcad file 121C3-02, read the introduction and move to page 2. 
Here a direction field has been set up for the function f(x,y) = e°** — 1; 
that is, for the differential equation 


er 
— —e 
a2 

Also, an initial condition is specified, as 


You studied this differential —]. 


equation in Activity 4.2(a). 


In the file, the initial values of | y = 0 when c= 0; that is, 2O7=6 


d denoted 
ee Note, to the right of the graph, that the value f(z, y,) = 1.718 is 
respectively by xp and yo. —. 
displayed, where n = 0. This is (to 3 d.p.) the value of 
{0,0 =e"? —1=e-1, 


which is the slope of the direction field at the starting point 

(Zo, yo) = (0,0). This starting point is denoted on the direction field by a 
small blue box, and the direction of the direction field at this point is 
coloured magenta. 


Ss 


You will now see how an approximate solution to this initial-value problem 
can be built up graphically, step by step. Ensure that all of page 2 from 
the heading ‘Solution curve’ to the bottom of the graph is visible on your 
screen. 


(a) Change the value under ‘Number of steps’ in turn to N = 1, N = 2, 
N =3, N=4 and N =5. In each case, observe the effects that the 
change in value causes. 


(b) Observe the effect of changing, in turn, the step size to h = 0.5 and 
the number of steps to N = 10. 


(c) Observe the effect of changing the initial values. For example, set 
Se = 1 ABS so = Z 


Comment 


© When ‘Number of steps’ is changed to N = 1, the first segment of the 
approximate solution curve is drawn, from the starting point 
(Xo, Yo) = (0,0) to the point (#1, y,) = (1, 1.718). The coordinates of 
these points appear in the tables to the right of the graph. The slopes 
of the direction field at these points, f(r, yo) = 1.718 and 
f (x1, y1) = 0.717, are also shown. 


The two broken blue lines that appear on the graph are temporary 
construction lines. They illustrate how the first segment of the 
solution curve is drawn as the hypotenuse of a right-angled triangle. 
The base (run) of this triangle is equal to the step size, which is 
specified above the graph as h = 1. The triangle is constructed, as 
shown in Figure 4.2, so that the gradient of its hypotenuse 1s equal to 
the gradient given by the direction field at its left-hand vertex. This 
gradient is f(xo, yo). On the other hand, the gradient (slope) of the 
hypotenuse is given as usual by the rule ‘slope equals rise over run’, 
where the run is the step size, h. The rise, which is the height of the 
triangle, is therefore given by 


rise = run Xx slope = hf (2o, yo). 


The coordinates (21, y,) of the top vertex of the triangle can then be 
calculated from the coordinates (Zo, yo), the rise and the run. This 
gives 


(21,41) = (o + run, Yo + rise) 


7 (x9 +h, Yo + hf (xo, yo)) 
-- (0 +1,0+ 1.718) = (1, 1.718). 


© When the number of steps is changed to N = 2, the second segment of 
the approximate solution curve is added, from the point 
(21, y1) = (1, 1.718) to the point (x2, y2) = (2, 2.4385). The broken blue 
lines now illustrate the ‘construction triangle’ for this second line 
segment, as shown also in Figure 4.3. The base (run) of this triangle is 
again equal to the step size, h = 1, but the gradient of its hypotenuse 
is equal to the gradient given by the direction field at (x,,y,). This 
gradient is f(x1, y,) = 0.717, as indicated by the table on the right of 
the screen. The coordinates (x2, y2) are therefore given by 


(v2, yo) = (a1 + run, y; +rise) 


= (4, th, y. +hf(x1,y1)) 
= (141, 1.718 + 0.717) = (2, 2.435). 


SECTION 4 Differential equations with the computer 


Comment on part (a) 
for NV = 1 


All of the values are shown to 
three decimal places. The 
slope of the direction field at 
(x1, y1) is now shown by a 
line segment on the graph, 
coloured magenta. 


(x1, 41) 


hf (xo, yo) 


Figure 4.2 First segment 


Comment on part (a) 
ma N= 2 


Figure 4.3 Second segment 


CHAPTER C3 


Comment on part (a) © 
lor N= 3, N = 4 and N =5 


Note that the values of 

f(%n, Yn) given in the table 
are negative for 

n = 2,3 and 4. These 
correspond to the negative 
slopes of the direction field at 
the corresponding points 

(Zn, Yn), as is clear from the 
graph. 


(2n41, Yn+1) 


Figure 4.4 (n+1)th 

segment ! 

Comment on part (b) © 
Comment on part (c) © 


With two steps now completed, it should be fairly clear how the 
process continues when the number of steps is changed to N = 3, 

N =4and N =5. The chain of line segments (which form the 
approximate solution curve) and the rows of corresponding values in 
the tables to the right of the graph build up one by one as N is 
incremented. Each time, the latest line segment is joined to the 
previous one, and its slope matches that of the direction field at its 
left-hand end. . 


Mathematically, this means that in constructing the line segment from 
c= Yn) to (2n41, Yn+1)s we have 
(Peat: Yn+1) ss (Se 5 2 run, Yn i rise) 


—_ (2. +h,yn + fiz, Yn)) 


This is illustrated in Figure 4.4. In other words, the sequences z,, and 
Yn are determined by the pair of recurrence relations 


(4 = 0) 1,2,=. 5); 


ba ae 4 25 a 


Yn4+1 = Un + hf (@n, Yn) 


as shown to the right of the graph on screen, together with the 
specified values for 2) and yo. These recurrence relations encapsulate 
Euler’s method. 


Lnt1 = In + h, 


Each value y,, is an estimate of the ‘true solution’ y at x = Z,; that is, 
Yn is an estimate of y(z,,). Clearly, the sequence of estimates obtained 
depends on the choice of both the step size, h, and the number of 
steps, NV. 


With WN steps of size h, Euler’s method provides NV + 1 solution 
estimates, spaced at regular horizontal intervals between the chosen 
starting value, x9, and 47) + Nh. 


Note that, when starting at (29, yo) = (0,0), either 5 steps of size 1 

or 10 steps of size 0.5 provide a final solution estimate at x = 5. By 
halving the step size and doubling the number N of steps, their 
product Nh, and hence also x) + Nh, remains the same. Notice, 
however, that the solution estimates obtained at x = 5 are different: 
the final coordinates are (25, ys) = (5, 0.986) in the first case and 

(710, Yio) = (5, 0.558) in the second. Halving the step size and 
doubling the number of steps improves the accuracy of the solution 
estimate. In the next activity you can investigate further how accuracy 
improves with reductions in step size. 


Setting 2» = 1 and yo = 2 produces a different chain of line segments, 
which is again an approximate solution curve for the differential 
equation, now with the initial condition y(1) = 2. (The same is true 
for other choices of initial values.) If h and N are unaltered from 
part (b), then the last row of values at the right of the screen 
corresponds to a point which lies outside the displayed graph region. 


Activity 4.3 showed how Euler’s method works. You saw a connected chain 
of line segments built up, one by one, as an approximation to a solution 
curve. However, this was a fairly inaccurate approximation. The 
construction utilises slope values provided by the direction field only at 
points which are a horizontal distance h apart, where the values assigned 
to the step size h in Activity 4.3 were first 1 and then 0.5. As a result, the 
approximate solution curve was based on very limited information. 
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More information can be extracted from the direction field by reducing the 
step size, provided that also the number N of steps is increased to 
maintain coverage of the x-values over which a solution is sought. ‘The next 
activity shows that such use of extra information leads to improvements in 
accuracy, and that an estimate for the solution y(x) at a particular chosen 
value of x can, in principle, be found to whatever accuracy is required. 


Activity 4.4 Using Euler’s method 


Move to page 3 of the document, and read the first paragraph. Then scroll 
down until all of the page from the heading ‘Graphical and numerical 
solution’ to the bottom of the graph is visible on your screen. 


The page is set up to solve the same initial-value problem that was 
considered in Activity 4.3(a) and (b), namely, 


dy 
— Se ae] O}-=5.4), 


The computation, using Euler’s method, takes place from the starting 
value x = 2p to the finishing value x = xval, where xval can be specified 
by the user. The step size h can be chosen as before. However, in contrast 
to the situation in Activity 4.3, it is not possible here to vary 
independently the number of steps, N. Given a step size h, the number of 
steps to be used is calculated automatically from the condition that the 
final step is to reach x = xval. As a result, the value calculated for yy is 
always an estimate for the value y(xval) of the true solution at x = xval. 


(a) Change the value of rval to 9 (from its starting value of 5). What 
changes on the screen, and what stays the same? 


(b) We now seek an estimate for the value of y(9) (the true solution value 
at z = 9). This is obtained by progressively reducing the step size. 


First note the value of yy that appears on screen (this is the estimate 
for y(9) obtained with step size h = 1, and hence with 9 steps). 


Now change the step size in turn to 0.5, 0.2, 0.1, 0.01, 0.001 
and 0.0001. In each case, note the corresponding value for yy. Use 
these values to estimate the value of y(9) to one decimal place. 


Solutions are given on page 34. 


Comment 


© As the step size h is decreased, the number of steps increases. With 
h = 0.0001, there are 90000 steps, for which the calculation may take 
an appreciable time on your computer. 


© With h = 0.5, the approximate solution curve still looks like a 
connected chain of line segments. However, with h = 0.2 the graph 
appears significantly more like a smooth curve, and this remains the 
case for smaller values of h. (The graph is still in fact made up of 
short line segments, but so is every graph drawn by Mathcad with the 
trace type set to ‘lines’ !) 


© Just as the numerical estimates for y(xval) appear to converge as the 
step size is decreased, so too do the approximate solution curves seem 
to converge towards a ‘limiting curve’ on the direction field. The 
graphs obtained provide increasingly close representations of the 
actual solution curve. 


You should still be working 
with Mathcad file 121C3-02. 


In other words, if h divides 
exactly into xval — xo, then 
N = (aval — zo) /h. 
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Euler’s method, as illustrated in Activities 4.3 and 4.4, was used to plot 
the solution curves in file 121C3-01, used for Activities 4.1 and 4.2. There, 
the number of steps was set to N = 1000, and calculations were made to 
draw the line segments to the left, as well as to the right, of the initial 
point. In that way, the line segments gave the appearance of a smooth 
solution curve over the whole horizontal graph range. 


If you are interested and have the time, then you might like to try the 
following optional activity, which involves a function f(z, y) that depends 
on the dependent variable y as well as on the independent variable x. For 
this particular example, it is also possible to check the outcome against a 
formula for the solution. 


Activity 4.5 Another initial-value problem (Optional) 


You should still be working (a) Use Euler’s method to estimate to three decimal places the value 

with page 3 of Mathcad of y(6), where y(x) satisfies the initial-value problem 

file 121C3-02. Before doing d 

anything else, reset the step = =7=—7—S, y(1) = 1. 

size to h = 1 (see the first dx 

note in the Comment below). (For the step size h, use in turn the values 1, 0.5, 0.1, 0.01, 0.001 
and 0.0001.) 


(b) Check that the function y = x — 4+ 4e*® is the solution to the 
initial-value problem in part (a). Hence find y(6) exactly. Does this 
agree with the value that you obtained using Euler’s method in 
part (a)? 


Solutions are given on page 34. 


Comment 


© If you do not start by setting h = 1, then any other change will cause 
recalculation for the most recently-used value of the step size, 


See the Comment for h = 0.0001, which may be time-consuming. Alternatively, recall that 
Activity 4.4 of Chapter B1, in any Mathcad calculation can be interrupted by pressing [Esc] and 
Computer Book B. then clicking ‘OK’ in the resulting option box. You may prefer to 


change here to ‘manual calculation mode’. 


© There is no need to alter any part of the file before page 3, nor to 
change the parameters which define the grid for the direction field. On 
page 3 you need to alter the direction field function definition to 
f(x,y) =x —y-—=3. Also, the values of x) and yo should both be set 
to 1, and the value of xval to 6. 


© Using the suggested values of fh in turn, the numerical estimates yy 
appear to converge and the approximate solution curves do likewise. 
The ‘chain of line segments’ is visible for h = 1 and for h = 0.5, but no 
departures from smoothness are apparent on the graph for smaller step 
S1Zes. 


Now close Mathcad file 121C3-02. 


‘Salsiens to (ote 


Chapter C1 


Solution 5.1 


The derivative of x?/5 is —0.4 when x = —1, and 0.2 
when z = 0.5. 


The function f(x) = x° — x* — 3x + 4 has derivative 
equal to zero at approximately z = —0.72 
and & = 1.39. 


Solution 5.2 


Where more than one expression is given below for a 
solution, the first is similar to the Mathcad output 
and the second is a form that you are more likely to 
obtain by hand. (You found each of the derivatives 
in parts (b)—(h) by hand in the main text of 
Chapter Cl, as indicated by the references below.) 
(a) 327 -—2x-—3 
(b) 4zr? 

See Activity 3.4(a). 


(c) (cos(t) — 2¢) _ (sin(t) — ¢*) 
exp(t) exp(t) 
cost — sint + t? — 2t 
= 
See Exercise 4.2(b). 
(d) 2cos(x”)z = 22 cos(x”) 
See pages 48-49. 
(e) —4sin(4z) 
See Activity 4.7(a). 
+3 
(f) 2tIn(t) + ait 
See Activity 4.2(c). 
1 In(u) u* +3 —2u? Inu 
(8) a ee 
u(u2 + 3) (u2 + 3) u(u? + 3) 
See Activity 4.4(b) 
exp(t) +1  e*+1 
exp(t) +t  et+t 
See Exercise 4.3(d). 


Solution 5.3 
Step 1: The stationary points of the function 


f(x) = 0.125\/1 + x? + 0.0625(2 — 2) 
are at x = £0.577 (to 3 d.p.). 


Step 2: Only the positive stationary point, 0.577, lies 
inside the interval [0,2]. The values of f at the 
interval endpoints are 


f(0) =0.250 and f(2) = 0.280, 


while the value of f at the stationary point in the 
interval is 


f (0.577) = 0.233 (all to 3 d.p.). 


Step 3: Hence the minimum value of f(x) for x in the 
interval [0, 2] is 0.233 at « = 0.577 (both to 3 d.p.). 


Hence the solution to the orienteer’s problem is to 
aim to join the path at a distance 0.577 km from the 
fixed point O on the path. 


(This agrees with the answer 0.58 km found in 
Activity 5.7(a) of Chapter A3 in Computer Book A.) 


Solution 5.4 
(c) Step 1: The stationary points of the function 
Vv 
sh 5+v+0.02v?2 


are at v = £15.811 (to 3 d.p.). 


Step 2: Only the positive stationary point, 
15.811, lies inside the interval [0,35]. The values 
of f at the interval endpoints are 


f(0)=0 and f(35) = 0.543, 


while the value of f at the stationary point in 
the interval is 


f (15.811) = 0.613 (all to 3 d.p.). 


Step 3: Hence the maximum value of f(v) for v 
in the interval [0, 35] is 0.613 at v = 15.811 
(both to 3 d.p.). 


(d) According to the model, a maximum traffic flow 
rate of about 0.61 vehicles per second can be 
achieved, by keeping the speed of traffic at 
about 16 ms! (that is, about 57 km per hour 
or 35 mph). 


= 


_ COMPUTER BOOK C 


Solution 5.5 


The answers are presented in a form as close as 
possible to the Mathcad output (which by default is 
given to 3 decimal places). Each answer agrees with 
that obtained by hand in the main text. 


Chapter C2 


Solution 5.2 


In each case, c is an arbitrary constant which has 


been added to the expression given by Mathcad. The 
answers to parts (a) and (b) agree with those 
obtained earlier by hand. 


1 
(a) / € + sa dx = In(x) + $ exp(3x) +c 
See Activity 1.2(a). 


3 1 
(b) / (= + ssin(6y) dy = one cos(5y) + c 
See Exercise 1.1(b). 


(c) [ic + cos(ax)) dx = ax + 


sin(ax) ie 2 


Solution 5.3 


In each case, c is an arbitrary constant which has 
been added to the expression given by Mathcad. The 
Mathcad answers in parts (a), (b) and (d) resemble 
closely those obtained by hand in the main text. 


(a) [(e-3)(@-1)ar = ig’ — 22° + 3r+ce 
See Example 2.1(a). 
or 
=" dx = 4259/2 —6/z+c 


Je 


See Example 2.1(c). 


(b) 


(c) [sw x dx = —4cos(x)sin(x) + $4+c 
In Activity 2.3(a), the answer $x — 4 sin(2x) +c 
was obtained by hand. The equivalence of this 
and the Mathcad expression follows from the 
double-angle formula sin(2z) = 2sin x cos z. 


tL 2 
See Activity 2.5(b) (iii). 


(e) [ue du = Zuexp(3u) — § exp(3u) +c 
(f) [@ In(5x) dx = $2" In(5zx) — 2° +c 


1 
(g) / ae dt = arcsin($t) +¢ 
The function arcsin is represented in Mathcad 
( 
by asin.) 


a2 


2 
(a) / e' dt = exp(2) — 1 = 6.389 
0 
See Activity 4.4(b). 
m/A 
(b) / (cos(5x) + 2sin(5zx)) dx = 4V2+2 


= 0.541 
See Exercise 4.1(a). 


2 
6 


See Exercise 4.1(b). (This answer is the 
Mathcad outcome from symbolic evaluation. ) 


(d) / e’ sint dt = 5 exp(m) + 5 = 12.070 
0 
See Exercise 4.1(c). 


Solution 5.7 


The answers are given by Mathcad to 3 decimal 
places. 


1 
(a) / a= 0.333 
0 


1 
(b) / V1— 2? dx = 1.571 
= 


SOLUTIONS TO ACTIVITIES 


Solution 4.2 


(a) (ii) The slope of the direction field at a point 
(x0, yo) appears to depend on the choice of x 
alone and not on that of yo. Correspondingly, 
there will be solution curves of just one type, 
with any two such curves differing only by a 
vertical translation. (This is a consequence of 
the fact that the function f(z, y) in this case 
depends only on 2.) 


Chapter C3 


Solution 4.1 


(a) The solution curves are described and sketched 
in the following table. 


(0,0) | Rough U-shape; 


steeper to right of 
minimum than to left. 
(0, 1) 
(0, 2) 


(b) The solution curves show three distinct types of 
behaviour: 


(iii) All choices for zo and yo give a solution 
curve which increases and decreases alternately, 
but has a rising trend. Each such curve may be 
obtained from that shown below by a vertical 
translation. ° 


Similar, but with 
minimum to the left 
and higher. 


Similar, but with 
minimum still further 
to the left and higher. 


Straight line 


psoas 1. Figure $3.1 


(b) (ii) The effect of the direction field on solution 
curves in this case is not so clear-cut. However, 
it appears that there may be three types of 
solution curve, as detailed below. 


(iii) There are solution curves that cross the 
x-axis and are increasing, as in the graph below. 
If zo = 0, then such curves are obtained by a 
choice of yo such that —1.25 < yo < 1.25. 


Downward curve with 
gradient decreasing. 


(i) Any solution curve through a point above 
the line y = —1 — x remains above that line. It 
decreases to a minimum and then increases. 


(ii) The line y = —1 — z is itself the solution 
curve through any point on that line. 


(iii) Any solution curve through a point below 


the line y = —1 — x remains below that line and 
decreases. 


The solution curves through the points (—3, —1), 
(—1,0) and (4,2) are of types (iii), (ii) and (i), 
respectively. 


Figure S3.2 


There are solution curves that lie above the 
x-axis, each having a minumum in the second 
quadrant, as in the graph in Figure $3.4. If 
Xo = 0, then such curves are obtained by a 
choice of yo such that yo > 1.26. 


a3 


Figure $3.3 


There are solution curves that lie below the 
x-axis, each having a maximum in the fourth 
quadrant, as in the graph below. If x» = 0, then 
such curves are obtained by a choice of yo such 
that yo < —1.26. 


Figure S3.4 


Solution 4.4 


(a) 


The graph extends to x = 9, but the portion to 
the left of s = 5 remains unchanged. In 
particular, the step size h does not alter, though 
the number JN of steps does (so that the graph 
reaches x = 9 rather than x = 5). 


A table of the values obtained is below. 


h 1 0.5 ° 0.22 -0.1 0.01- 0.00] GOGH 


yn 3.916 3.347 3.002 2.887 2.783 2.772 2.771 


The values of yy appear to be converging and, 
given the level of agreement between the last 
two estimates, it seems reasonable to deduce 
from them an estimate for y(9) that is accurate 
to one decimal place, that is, y(9) = 2.8. (In 
fact, it looks likely that y(9) = 2.77 to 2 decimal 
places.) 


Solution 4.5 


(a) 


The table below gives the values obtained using 
Euler’s method, from file 121C3-02. 


h 1 0.5 O.1 0.01 0.001 0.0001 
yn 2.000 2.004 2.021 2.026 2.027 2.027 


The agreement of the last two estimates suggests 


that y(6) = 2.027 to three decimal places. 


For the given function y = x — 4+ 4e!~*, we 
have 


y(1) =1—4+4e° =1, 


so that the initial condition of the problem in 
part (a) is satisfied. 


The derivative of the given function is 
dy _ 
dx 

whereas we have 
e—-y—-3=2-—(r—4+4e' 7) -—3 


=1-4e!-*, 


ee ee 


Hence the given function also satisfies the 
differential equation from part (a), and therefore 
satisfies all the conditions of the initial-value 
problem. 


The value of this solution function at x = 6 is 
y(6) =6—4+4+4e' °=2+44e > ~ 2.027. 


This value agrees to 3 decimal places with that 
obtained using Euler’s method. 
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